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We present numerical evidence that the SU(4) pure-gauge dynamics has a finite-temperature first-order phase 
transition. For a 6 x 20^ lattice, this transition occurs at the inverse-square coupling of S/g^ ~ 10.79. Below this 
and above the known bulk phase transition at 8/(?^ ~ 10.2 is a confined phase in which we find two different string 
tensions, one between the fundamental 4 and 4* representations and the other between the self-dual diquark 6 
representations. The ratio of these two is about 1.5. The correlation in the adjoint representation suggests no 
string forms between adjoint charges. 
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There are renewed interests in SU(A'c) pure 
Yang-Mills theory with large Nc: 

1) Finite-temperature phase structure of quan- 
tum chromodynamics (QCD) would be easier to 
understand if the SU(Afc) pure Yang-Mills system 
has a second order phase transition for iVc > 4 
|l] . With standard large- iVc analysis where Ncg'^ 
is held fixed, the 'Zj{Nc) deconfinement transition 
occurs at Td ^ 0(1), separating confining phase 
with free energy F ^ 0(1) and dcconfining phase 
with F ~ 0{N^). The deconfining temperature 
T<j 0(1) is not efi'ected if TV/ and g^N^ are held 
fixed and Nc — > oo. If the transition is first or- 
der, it is not effected either. So large N^, is not a 
reasonable guide for T ^ QCD phase structure 
with Columbia phase diagram 1^, unless SU(Afc) 
pure Yang-Mills dynamics has second order de- 
confining phase transition for all Nc > 4. 

2) New developments in M/string theory 
predict such things as glueball spectrum at large 
Nc and large or ratio between different string 
tensions for A^c > 4 [Q . 

3) The dimensionless ratio T^j \fo of the de- 
confining temperature Td and string tension cr 
is expected to be independent of Nc with a 
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value Y/377r(I?^2) with D being the space-time 
dimensions |^]. 

Here we report the results of our numerical 
investigation on the order of deconfining phase 
transition and the ratio of string tensions for 
Nc = 4 dj. We use single-plaquette action de- 
fined in the fundamental 4-representation of the 
SU(4) gauge group. Combinations of pseudo- 
heatbath or Metropolis and over-relaxation al- 
gorithms are used in updating 4, 6 and 8 x8'^, 
12^, 16'^ or 20'^ lattices. Various worksta- 
tions are used for the numerical calculations, 
while migration to the RIKEN BNL QCDSP 
mother boards is planned. We look at the fol- 
lowing observables: plaquette, Polyakov loops, 
L[x) = {l/Nc)trllf^^U{x,t;i), in 4 (fundamen- 
tal), 6 (antisymmetric diquark), 10 (symmet- 
ric diquark) and 15 (adjoint) representations, 
deconfinement fraction, and Polyakov loop cor- 
relation (L(0)L(f)*) - r-iexp(-F(r)/r) - 
exp(— Lfcrr — Inr) in 4, 6, 10 and 15 representa- 
tions. 

This SU(4) pure Yang-Mills system is known 
to have a bulk phase transition near f3 = 8/g'^ ~ 
10.2, separating two confining phases 0: across 
this transition the plaquette jumps discontinu- 
ously but the average Polyakov line in the fun- 
damental 4 representation remains zero on both 
sides. However if the lattice extent in tempera- 
ture direction Lt is too small this bulk transition 
drives a first-order finite-temperature deconfining 
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Figure 2. Polyakov line correlation in 4 (+) and 
6 (x) representations on a 6 x IG'^ lattice at /3 = 
10.70. From the slopes we find different string 
tensions for these representations. No signal was 
obtained for 10 and 15 representations, probably 
because the coupling is too strong. 



Figure 1. Time histories of the fundamental 
Polyakov loop magnitude and argument (in units 
of tt) from a 6 X 20^ lattice at /3 = 10.79 (above) 
and magnitude histogram (below). Confined and 
deconfined phases coexist at this temperature 
suggesting a first-order deconfining phase tran- 
sition. 



transition just above itself in /? We confirmed 
this is indeed the case for Lt = 4. But for Lt > 6 
we have good enough separation between the bulk 
and finite-temperature transitions as desired. 

As is shown in Figure |l], on a 6 x 20'^ lat- 
tice we confirmed coexistence of confined and de- 
confined phases at temperature /3=10.79: This 
strongly suggests a first-order deconfining phase 
transition. Work in progress confirms this phase 
coexistence as we extend the simulation from the 
present 3500 evolution steps (1 evolution — 5 heat 
bath 4- 1 over relaxation steps) to 20000 steps. 
We plan further study with finite-size scaling. 



String tensions in SU(iVc) pure Yang-Mills sys- 
tem is classified by its center Z{Nc) iVc-ality. 
With Nc — 4, the fundamental (4) charge has 
4-ality k = 1, the two diquark (6 and 10) charges 
k = 2, and adjoint (15) fc = 0. The string ten- 
sions between these charges and their anticharges 
are predicted to behave as ak oc minjfc, A^c — k} 
by a standard strong-coupling analysis, k{Nc — 
k) by another strong coupling analysis Q, and 
sin(fc7r/Ac) by a SUSY strong coupling analysis 
0. Generally the ratio cTfc/tTi should fall in the 
interval 1 < fifc/cri < 2 §]. Note also that TVc = 4 
is the first example with different string tensions: 
in SU(3) pure Yang-Mills system the fundamen- 
tal (3) and the symmetric diquark (6) tensions 
are the same ||lO| . 

In our numerical calculation on a 6 x 16'^ lat- 
tice at /3 = 10.70 (see Figure ||): we find a clear 
difference between 4- and 6-string tensions ex- 
tracted from 4- and 6-Polyakov loop correlations. 
From fitting these data we have CT4=0. 068(4) and 
(T6=0.108(17). At a stronger couphng of (3=10.65 
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Figure 3. Polyakov line correlation in 4 (+) , 6 ( x ) 
and 15 (*) representations on a 8 x 12^ lattice at 
/3= 10.85. The adjoint (15) signals now suggest 
there is no string for this representation. From 
the slopes of the former two lower-dimensional 
representations we confirm different string ten- 
sions for them, and by comparison there is no 
tension seen in the adjoint representation. 



their values are 0.086(3) and 0.142(57) respec- 
tively. Thus their ratio ae / does not show much 
temperature dependence and falls in the interval 
(1,2) as it should. We are yet to see any signal 
for 10 and 15 representations from this lattice, 
probably because of too strong coupUngs. On the 
other hand at a weaker coupling of /3=10.85 on a 
smaller lattice of 8 x 12^ (see Figure |^) we find 
clear flattening of the adjoint (15) correlation. 

For thermodynamics of SU(4) pure-gauge the- 
ory we confirmed that the bulk transition and 
T ^ phase change are separated on Lt > 6 
lattices, the bulk transition is at f3i, ^ 10.2 and 
the finite-temperature Ist-order phase transition 
is at /3d ~ 10.79 (L* = 6) and > 10.9 (Lt = 8), 
and easier to establish than weakly first-order 
SU(3). For string tensions, at Lt=6 we find sig- 
nals for different string tensions in fundamental 
4 and antisymmetric 6 representations, but no 
signal yet for symmetric diquark 10 and adjoint 



15 representations. These strings satisfy a rela- 
tion, 1 < (J(i/a4 < 2, as they should, and the 
ratio does not show any strong temperature de- 
pendence. Combining these findings for thermo- 
dynamics and string tensions at Lt=6, we find an 
inequaHty: Td/^a^iT^ O) < Tdl^a^i^T ^ Td) 
- 0.64 < ^3/7r(i:)-2), just like in SU(2) and 
SU(3) pure-gauge results. At a weaker coupling 
oi ji — 10.85 using a Lf=8 lattice we now have 
rather good evidence that there exists no string in 
the adjoint (15) representation. We plan further 
investigation on larger and finer lattices, probably 
using smaller partitions of the QCDSP parallel 
supercomputer at RIKEN BNL Research Center. 
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